Generalization of escape rate from a metastable state driven by external 

cross-correlated noise processes 
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We propose generalization of escape rate from a metastable state for externally driven correlated 
noise processes in one dimension. In addition to the internal non-Markovian thermal fluctuations, 
the external correlated noise processes we consider are Gaussian, stationary in nature and are of 
Ornstein-Uhlenbeck type. Based on a Fokker-Planck description of the effective noise processes with 
finite memory we derive the generalized escape rate from a metastable state in the moderate to large 
damping limit and investigate the effect of degree of correlation on the resulting rate. Comparison 
of the theoretical expression with numerical simulation gives a satisfactory agreement and shows 
that by increasing the degree of external noise correlation one can enhance the escape rate through 
the dressed effective noise strength. 
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I. INTRODUCTION 

The theory of fluctuation induced barrier crossing 
dynamics was first discussed in the seminal paper by 
Kramers [l[ where he considered a thermalized Brownian 
particle trapped in a one dimensional well separated by 
a barrier of finite height from a deeper well. The particle 
was supposed to interact with the environmental degrees 
of freedom that act as a thermal reservoir. The thermal 
environment exerts a damping force on the particle but 
simultaneously thermally activates it so that the parti- 
cle effectively gains enough energy to cross the barrier 
of finite height. Since its inception, Kramers' model has 
been extensively revisited both theoretically [1, d, 0| and 
experimentally 0, @, 0] constituting a vast body of liter- 
ature Hi]. 

One of the essential conditions of the traditional ap- 
proach to study activated rate processes is the main- 
tenance of a balance between the two opposing forces, 
the thermal fluctuations and the dissipation, which the 
Brownian particle experiences while in contact with the 
thermal bath. Since these two counter balancing forces 
have a common origin, the heat bath, it can be shown 
easily that they are connected by fluctuation-dissipation 
relation (FDR) [ToJ ■ A typical signature of FDR is that 
in the long time limit the Brownian particle attains an 
equilibrium Boltzmann distribution for an initial canon- 
ical distribution of the bath's degrees of freedom [ll[ • In 
nonequilibrium statistical mechanical terminology, such 
systems are referred to as closed system [l2j]. It may 
happen sometime that an additional source of energy in 
the form of fluctuations can be pumped from outside for 
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which there is no counter balancing force like dissipation 
EI El EH- In absence of FDR, such external in- 
put of energy makes the system open and in contrast to 
the closed system the equilibrium Boltzmann distribu- 
tion gets replaced by a steady state distribution (SSD) 
in the long time limit. It may therefore be anticipated 
that the absence of FDR tends to make the SSD func- 
tion dependent on the strength and correlation time of 
external noise as well as on the dissipation of the system 
EH EH • It is pertinent to point out that though ther- 
modynamically closed systems with homogeneous bound- 
ary conditions possess, in general, time-dependent solu- 
tion, the driven open systems may settle down to compli- 
cated multiple steady states when one takes into account 
nonlinearity of the systems. 

The origin of the noise in the thermodynamically open 
system driven by two or more random forces may be dif- 
ferent. The barrier crossing dynamics with multiplica- 
tive and additive noises arose strong interest in the early 
eighties. The noise forces which appeared in the dynam- 
ical system were usually treated as random variables un- 
corrected with each other. However, there are situations 
where noises in some open systems may have common 
origin. If this happens then the statistical properties of 
the noises should not be very much different and can 
be correlated to each other. The cross-correlated noises 
were first considered by Fedchenia [I?) in the context of 
hydrodynamics of vortex flow. The interference of addi- 
tive and multiplicative white noises in the kinetics of the 
bistable systems was analyzed by Fulinski and Telejko 
fl8l |. Maduriera et al. [l9| have pointed out the prob- 
ability of cross-correlated noise in ballast resistor model 
showing bistable behavior of the system. Mei et al. [2(| 
have studied the effects of correlations between additive 
and multiplicative noise on the relaxation of the bi-stable 
system driven by cross-correlated noises. It is now well 
accepted that the effect of correlation between additive 
and multiplicative noise is indispensable in explaining 
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phenomena like phase transition, transport of motor pro- 
teins, etc [21] . As the presence of the cross-correlated 
noises changes the dynamics of the system l 22|, it is 
expected that there may exist some additional effect of 
cross-correlation on the barrier-crossing dynamics. The 
study of this additional effect acted precisely as the cat- 
alyst that triggered the present study. From our formal 
development and corresponding numerical application, it 
will be revealed that the strength of correlation of noise 
process has pronounced effect on the behavior of the es- 
cape rate of barrier crossing process. 

We are now in set to discuss the physical motivation 
of our model. One can think about a system, simulta- 
neously coupled to two different heat baths. The two 
heat baths are in turn driven externally by a random 
force e(t) (say). The response of e(t) to the two baths 
are different. Now, if the system-bath coupling be lin- 
ear, apart from thermal noises due to the presence of 
heat baths, the system will encounter two additive dressed 
noises (dressed due to the bath-external noise coupling). 
As both the baths are modulated by same noise with 
different response, the dressed noises will be mutually 
correlated. Another simple model can be invoked by 
considering the following system-reservoir Hamiltonian: 
H = H s + H SB + H int , where H s = p 2 /2 + V(x) is the 
system Hamiltonian and Hsb = Ylj [p] + w f ("ft ~ c j x ) 2 }/^ 
being the bath Hamiltonian (with interaction term be- 
tween the system and the heat bath), where the heat 
bath is assumed to be consisting of harmonic oscillators 
with characteristic frequency set {ujj} as well as {qj,Pj} 
are bath variables, and H lnii = £\ [Kjqje(t) + \jp -e(t)]. 
While writing the Hamiltonian we have considered unit 
mass of the system and the bath oscillators. In this 
model, initially the bath is in thermal equilibrium at the 
temperature T in the presence of the system, and then 
is externally modulated by a noise e(t). The coupling 
between the bath and the noise is such that it linearly 
excites both positions qj and momenta Pj of the bath 
with the different coupling constants Kj and Xj respec- 
tively. If we now construct the equation of motion for the 
system variable, it can be easily shown that the equation 
will contain, apart from the thermal noise, two additive 
mutually correlated noises. As a more physically moti- 
vated example, one can consider a typical photochemical 
reaction where both reactant and solvent are exposed to 
a weakly fluctuating light source. In such a case, the 
reactant will be driven by two noises, first the external 
noise due to the fluctuating light source and second, due 
to the presence of the solvent which is also driven by the 
light source [l6j]. These two noises will also be corre- 
lated as both the reactant and the solvent are exposed 
to the same fluctuating light source. The above men- 
tioned physical examples lend firm support to our study 
presented in this article. 

The organization of the present work is as follows. 
In the next section (Sec. |TTJ) we briefly describe the 
phenomenological model for a system driven by cross- 
correlated external fluctuations and characterize the sta- 



tistical properties of external driving. In Sec. IIIII we de- 
rive the generalized escape rate from a metastable state 
and show different limiting cases in moderate to large dis- 
sipation regime. Computational details and results are 
provided in Sec lIVI The paper is concluded in Sec. [V] 

Before embarking on a discussion of the development 
and application of our theory for escape rate from a 
metastable state driven by external cross-correlated noise 
processes, in the next section, we will discuss the essen- 
tial ideas of our model. This will motivate us towards 
the types of physically appealing approximations needed 
in generating the required rate equations. 

II. THE MODEL 

To start with we consider the motion of a Brownian 
particle of unit mass moving in an external force field 
V(x). In the course of its dynamics the particle experi- 
ences a random force f(t) as well as a counter balancing 
frictional force "/(t) both originating from the immediate 
thermal environment, the heat bath, to which the parti- 
cle is in contact. Since the aforesaid forces have a com- 
mon origin they are connected by the FDR (see Eq.(2) 
below) flOI | . Apart from the internal thermal noise, we 
assume that the particle is driven by two external non- 
thermal, stationary, Gaussian, Ornstein-Uhlenbeck noise 
processes e(t) and ir(t) both of which are correlated to 
each other by the correlation parameter A. The dynamics 
of the particle can then be described by the generalized 
Langevin equation 

x + [ ~f(t-t')x(t')dt' + V'(x)=f(t)+e(t)+Tr(t). (1) 
Jo 

Here the friction kernel j(t) is connected to the internal 
noise f(t) by the FDR of the second kind [l(| 

(f(t)f(t'))=k B T 1 (t-t') (2) 

where T is the thermal equilibrium temperature and k B 
is the Boltzmann constant. The form of the Langevin 
equation (TT|) we have considered here guarantees that the 
nature of the dynamics is of non-Markovian type and it 
is due to the finite correlation effect of the thermal en- 
vironment on the Brownian particle. The nature of the 
dissipation kernel j(t) is very much dependent on the na- 
ture of the coupling of the particle to the heat bath and 
on the distribution of the bath modes [ll[ . In the limit 
of vanishing correlation effect and for instantaneous dis- 
sipation, i.e. for Markovian dynamics, Eq.(2) reduces to 
(/(*)/(*')) = 2>yk B T5(t-t') with 7 being the dissipation 
constant. 

At this point it is pertinent to mention that by the 
application of an external random force, fluctuations are 
created in a deterministic system. As example, one may 
cite a noise generator inserted into an electric circuit or a 
growth of species under the influence of random weather. 
For such cases, the external fluctuating force is never 



3 



completely (S-correlated or white. On the other hand, 
the internal or intrinsic noise is caused by the thermal 
fluctuations created due to the coupling with the envi- 
ronmental degrees of freedom and cannot be completely 
switched off. In general, the physical processes concern- 
ing chemical reactions, growth of population etc. are of 
later type. For an illuminating discussion in this con- 
text we refer to the book by van Kampen (23|. At this 
point we note that a system where internal noise is always 
present, can be driven by external noise also. From a mi- 
croscopic point of view, the above equations |T]) and (J2|) 
can be derived from a Hamiltonian, where the system 
is coupled with a heat bath consisting of a set of har- 
monic oscillators with different characteristic frequencies 
[ill ] and is simultaneously driven by two external, cross- 
correlated noises, e(t) and 7r(i). The system-bath inter- 
action generates the internal noise f(t), statistical prop- 
erties of which will depend on the frequency spectrum 
of the heat bath and on the nature of the system reser- 
voir coupling. For a finite correlation time r c , (which can 
be considered as inverse of the cut-off frequency of the 
bath), the internal noise will be colored, while for r c — > 0, 
the noise will be ^-correlated or white. Furthermore, the 
nonlinear system-reservoir coupling yields multiplicative 
noise while bilinear coupling gives additive noise. For 
further discussion, we refer to the book of Lindenberg 
and West [12j. In our present study, all the noises have 
been assumed to be colored and additive. 

The external noise processes we have considered are 
independent of the dissipation kernel, hence there exist 
no corresponding FDR for them. In addition to that, we 
also assume that the external noises do not influence the 
internal noise process and hence e(t) and 7r(£) are sta- 
tistically independent of f(t). The physical situation we 
address here is that the system is in thermal equilibrium 
at t = in presence of the thermal bath but in absence 
of the external noise processes. At t = 0+, the external 
fluctuations are switched on and the system is driven by 
the two external noises e(t) and n(t) which are correlated 
to each other. The system dynamics is then governed by 
the generalized Langevin equation ([1]). We characterize 
the statistical properties of the correlated external fluc- 
tuations by the following set of equations 
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<e(t)e(t')> 
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In the above equations, D e and r c are the strength and 
correlation time of the noise e(t), while D n and r n cor- 
respond to the noise n(t). A denotes the degree of corre- 
lation between the noise processes e(t) and 7r(t). At this 



point, we define [£(£) = (e(i) +ir(t)) ] as the effective 
external noise process whose statistical properties can be 
defined using the properties of e(t) and ir(t). Our follow- 
ing analysis will be based on this effective noise £ (t) and 
we shall study the dependence of physical quantities in 
terms of this effective noise £(£) and the effective corre- 
lation time. Since both the noise processes are Gaussian 
and simultaneously stationary with zero mean, £(i) will 
also be stationary Gaussian with zero mean. Further- 
more, the second moment of £(t) is given by 



TR V T R 



(4) 



where the strength Dr and the correlation time tr of the 
effective noise can be written as 

/>OC 

Dr = / di(£(*)£(0)) = D e + 2X^D € D n + D v (5) 
Jo 

D R Jo 

(6) 



D F 



D e r e + 2\y/D e D^r + D v Tn 



It is apparent from equations ((4j[6j) that variation of the 
degree of correlation A will change the strength and cor- 
relation time of the effective noise. 

From the above equations (|4][6]), it is clear that because 
of the cross correlation the Brownian particle realizes the 
effect of two external noise processes by the effective noise 
process £(t) with the effective noise strength Dr and the 
noise correlation time tr containing the noise strength 
and correlation time of e(t) and ir(t) and the degree of 
noise correlation parameter A. 



III. GENERALIZATION OF ESCAPE RATE 

The modifications of the standard escape rate from 
a metastable state in presence of correlated fluctuations 
are twofold. First, in presence of the external correlated 
fluctuations, the dynamics around the barrier top gets 
affected in such a way that the stationary flux across the 
barrier top gets modified. Second, the equilibrium Boltz- 
mann distribution at the source well gets replaced by a 
steady state distribution reflecting the signature of extra 
energy input due to the open system. With these major 
changes in the dynamics we then derive the generalized 
escape rate and show various limits of the rate expression 
in the following analysis. 

To analyze the dynamics across the barrier top we first 
linearize the potential V(x) around the barrier top at 
x w 



V(x) = E b 



1 



2 2 , 



> 



(7) 



where Wh is the linearized frequency at the barrier top and 
Eb (= V(0)) is the barrier height. Thus, the linearized 
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version of the Langevin equation (fTJ) takes the form 



From Eqs.(|l7[) and (fT9|) we observe that 



x + 7 (t - t')x{t')dt' - cu^x = F{t) (8) 
Jo 

where F(t) is defined as 

= /(*)+£(*). (9) 
The general solution of Eq.Q is given by 



x{t) = (x(t)) + / M b (t - t')F{t')dt (10) 
Jo 



where 



(x(t)) =v M b (t)+x xi(t) 



(11) 



with xq — x(0) and vq — x(0) being the initial position 
and velocity of the particle and 



X|(t) = l+^ / M b (t')dt'. 
Jo 

The kernel M b (f) is the Laplace inverse of 

1 



M b (s) 



s 2 + S7(s) - uj'I 



(12) 



(13) 



with 7(s) = J °° exp(— st)j(t)dt, being the Laplace trans- 
form of j(t). The time derivative of Eq. lfTO]) gives 



v (t) = («(*)) + / m 6 (t - t')F{t')dt' (14) 
Jo 



where 



and 



(v(i)) = v m b (t) + u%x M b (t) 
dM h (t) 



m b {t) 



tit 



(15) 



(16) 



As both the internal and external noise processes are sta- 
tionary, one can explicitly use the stationary property 
of these fluctuations to write the correlation function of 
F(t) as (F(t)F(t r )) = c(t - t') along with the symmetry 
condition c(t — t 1 ) = c(t' — t). Using this form of station- 
arity we calculate the variances in terms of M b (t) and 
m b (t) as 

°L(t) = ({x(t) - (x(t))} 2 ) 

= 2 I Ahit^dh I M b {t 2 )c{t 1 -t 2 )dtin) 



*l(t) = (Ht)-(v(t))} 2 ) 



= 2 / m 6 (ti)^i / m b (t 2 )c{tx - t 2 )dt 2 (18) 
Jo Jo 

<&,(*) = ([*(*) - (x(t))]Ht) - <»(*)>]> 

m 6 (ti)dti / M b {t 2 )c{t x -t 2 )dt 2 (19) 
o Jo 



5 (*)/2. 



(20) 



Using the method of characteristic function [24| we then 
write the Fokker-Planck equation for probability distri- 
bution function P{x, v, t) near the barrier top as 



dP dP 2/ , dP , MvP) 

± , ,d 2 P , . . d 2 P 

-Mt) — + Mt)- 



dv 2 



dxdv 



(21) 



where the subscript '6' signifies the dynamical quantities 
defined at the barrier top and they are given by 



%(t) = -j t ^Y b (t) 



-M b (t)m h (t) + M 2 (t) 
Y b (t) 

mb(t) + f M b {t')dt'\+M 2 (t) 



nw = - 



Mt) = u 2 b {t)al x {t) + %(t)a 2 xv (t) + & 2 xv (t) - a vv (t) 2 . 

(22) 

Although bounded, these time-dependent parameters 
may not always provide long time limits. In general, 
one has to work with frequency ui b {t) and friction %(t) 
for analytically tractable models. In absence of external 
noise <f> b (t) = k b T%(t) and ip b (t) = {k B T / u) b )[Q b - cjg] 
and in the Markovian limit % (t) = j5(t), Q b — w b and 
consequently, ip b (t) vanishes. 

To calculate the stationary distribution near the bot- 
tom of the left well, we now linearize the potential V(x) 
around x rj Xq. The corresponding Fokker-Planck equa- 
tion describing the dynamics at the source well can again 
be constructed using the method of characteristic func- 
tion, 



dP 
~dt 



with 



dP , dP , ,d(vP) , , fi 2 P 

-v— - wg(*)x-5- + 7oW^r^ + Mt) 



dv 



dv 



+Mt) 



d 2 p 

dxdv 



dv 2 
(23) 



70 (i) = ~J t ^Yo(t) 



M (t)m (t)+M 2 (t) 
Yo(t) 

mo{t) W + oj 2 f M Q {t')dt'\+M*[t) 



u, 2 (t) = 
Yo(t) = 



Mt) = u 2 (t)a 2 xx (t)+MtW 2 vv (t) + ^vvit) 

Mt) = ^l{t)a 2 xx (t)+Mt>L{t)+o 2 xv {t)~a 2 vv {t). 



(24) 
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Here, the subscript '0' signifies the dynamical quanti- 
ties corresponding to the bottom of the left well. As 
mentioned earlier, in absence of external noise and in 
the Markovian limit, ipo(t) vanishes with 7o(t) = jS(t), 
</>()(£) = iksT and ujq — luq. As a result, we recover the 
Fokker-Planck equation for a harmonic oscillator with 
frequency luq. Thus, we identify Eq. (|23|) as the general- 
ized version of the non-Markovian Fokker-Planck equa- 
tion for a harmonic oscillator which is driven by two 
externally correlated noise processes. At this juncture 
it is pertinent to mention the fact that the form of 
our Eq. (f2"3")) is exactly identical with that of the non- 
Markovian Fokker-Planck equation for harmonic oscilla- 
tor derived earlier by Adelman 25]. The steady state 
[dP(x w xq, v)/dt ~ 0] solution of Eq. ([23|) is given by 



p st(x,v) = ^exp 



LUqX 



2Dn 



Dn 



4>o 



(25) 



where Dq — 0o/7o- The quantities ipo, 4>o and 7o are long 
time (steady-state) limit of the time dependent functions 
ipo(t)i 4>a{t) and 70 (t) respectively and Z is the normal- 
ization constant. The above solution (|25[) can be veri- 
fied by directly putting it into the steady state version 
(dP/dt = 0) of Eq.(H21). The distribution function 
is the steady state counterpart of equilibrium Boltzmann 
distribution (exp{ — [v 2 +V (x)]/ kBT}) for nonequilibrium 
open system. In the limit of pure thermal processes, i.e. 
in absence of external fluctuating driving force, it is easy 
to recover the equilibrium Boltzmann distribution from 
([23)) . The justification for using the distribution (|25|) is 
the following. In the traditional theory of activated rate 
processes within the framework of pure thermal fluctua- 
tions the equilibrium Boltzmann distribution is necessary 
to initially thermalize the reactant state [J, Q ■ On the 
other hand, for nonequilibrium open system, a constant 
input of energy through external fluctuating driving force 
forbids the system to attain the equilibrium state, hence 
the system approaches (as in our case) towards the steady 
state (|23|). an analogue of equilibrium state, to initially 
energize the reactant state by the effective temperature 
like quantity (a complex function of 70, (po and ipo) em- 
bedded in the distribution function (f2"5)) . 

To calculate the stationary current across the barrier 
top (for pure thermal processes) within the framework of 
Kramers' original reasoning one considers an equilib- 
rium Boltzmann distribution multiplied by a propagator 
G(x,v) and uses it to solve the Fokker-Planck equation 
(|2ip . From the reasoning of the previous paragraph, it 
is clear that in our analysis, the equilibrium distribution 
should be replaced by a SSD that contains all the infor- 
mation about the potential around the barrier top and 
the effective temperature like quantity for nonequilibrium 
open system. In the limit of pure thermal processes, i.e. 
in absence of any external fluctuations the SSD should 
reduce to the equilibrium Boltzmann distribution. Thus, 
in the light of Kramers' ansatz we consider a solution of 



Eq. (f21j) at the stationary limit as 

V(x) 



P st( X i V ) = ex P 



V 

2D~b 



G(x,v) (26) 



where D b = (pb/% and ip b are the long time limits of 
the corresponding time-dependent quantities calculated 
at the barrier top region and V{x) is the linearized poten- 
tial near the barrier top (x « 0) with a renormalization 
in its frequency, 



V(x) ~E b - \u> 2 b x 



(27) 



where uj b is the long limit of Co b (t) and E b is the barrier 
height. In writing Eq. (|2"T)) it has been assumed that the 
position of the maxima of the potential V(x) and the 
barrier height remains unchanged while considering the 
memory effect in the dynamics. The non-Markovian ef- 
fects are reflected only in the frequency. The ansatz of 
the form ([26]) denoting the SSD is motivated by the local 
analysis near the source well and the top of the barrier 
in the Kramers' sense. Inserting Eq. (f2"6"l) in Eq. (|2~lj) we 
get in the steady state (dP/dt = 0) an equation for the 
function G(x, v) as 



1+^ 
8 2 G 



dG 

dx 



Db + tp, 



-uj b x + -y b v 



dG 
dv 



Wb „ „ = 0. 



dv 2 ' ^" dxdv 
We now introduce a variable u as 

u = v + ax 



(28) 



(29) 



where a is a constant to be determined. With the help 
of the transformation (|29l) , Eq. (|28|) reduces to 



(0 b + atpb) j. 2 



fG 
du 2 



D 



Db + tpb 



b -2 
UJ b X 



+ 1 76 + a ^1 + 
At this point, we define 



4>b 

D b 



dG 
du 



= 0. 



(30) 



Db + ip b 



uj b x + < 7(, + a I 1 + 



ipb 
D b 



-jiu (31) 



where fi is another constant to be determined. From 
Eqs. (f2"9")l and (|31[) . we find that the constant 'a' has two 
values, 



with 



a± 



ipb 



B 
2A 



± 



B 2 



C 



AA 2 A 



(32) 



A = ' 1 + if ' ' B = % and ° = 



(D b + rjft) 
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With the help of Eq.([3Tj), Eq.j30J) can then be written as 

(33) 



dw z au 



where 



A 



4>b + ai>b 
The general solution of Eq. (|33|) is 

1 



G(u) = G 2 I exp 
'o 



-Az 2 



dz + Gi 



(34) 



(35) 



where G\ and G 2 are the constants of integration. We 
look for a solution which vanishes for large x. This con- 
dition is satisfied if the integration in Eq.l[35|) remains fi- 
nite. It is easy to understand that the integral in Eq. (j3"5|) 
converges for \u\ — » oo if and only if A is positive. The 
positivity of A depends on the sign of 'a' and we observe 
that the negative value of 'a\ i.e., 'a_' guarantees the 
positive value of A. 

To determine the value of G\ and G2 , we now demand 
that G{x, v) — > for x — > +00 and for all v. This condi- 
tion yields G\ — G 2 y / n/2A, so that 

"/tt\ 1/2 r ( 1 



G(u) = G 2 



\2AJ 



exp 



Az< 



(36) 



Consequently, the stationary solution near the barrier top 
becomes 



i$(a:wO,t>) 



G2 exp 



D b + %p b 



(a) 



1/2 



exp 



2^ 



+ cjr(cc ~ 0, i>) exp 



2A 



(37) 



with 



ff(a;,u) 



rV+a-X /I \ 

y exp ( --Az 2 J dz. (38) 



Since the steady state current j across the barrier is de- 
fined as 



j = j vP b st {x w 0,v)dw, 



(39) 



we get using Eq. (f3"Tj) 



j = G 2 D b 



2n 



1/2 



t 1 exp , 

A + D7 V V Z? h + V'b 



£5 



(40) 



To obtain the remaining constant G 2 , we note that as 
x — > —00, the term G 2 [^Jtt/2A + g] in Eq. (|3"6")) reduces 



to G2 \/7r/2A. We then obtain the reduced distribution 
function 1261 in 2: as 



P s b t (x - -co) = 2^G 2 [ ) exp 



1/2 



£>0 + Vo 



. (41) 

Similarly, we derive the reduced distribution function in 
the left well around x ~ xq as 



1 



P° (2) = -y/2nD a exp 



E 



D + ipa 



(42) 



where we have used Eq. (|25p and employed the expansion 
of V(x) as 

V(x) ~ E + ]^l{x - x f , x w z . (43) 

At this point, we impose another condition that at 
x ?s xo, the reduced distribution function (|4ip must go 
over to the stationary reduced distribution function (|4"2|) 
at the bottom of the left well. This matching condi- 
tion [3, [III, [l6| along with the normalization condition, 
f_°° Pg t (x,v)dxdv — 1 gives the value of the remaining 
constant G 2 as 



1/2 



w 



TV [87T 3 (Z? + ^o)] 1/2 



(44) 



Hence, from Eq. (|4"0"]) . we get the expression for the nor- 
malized current or barrier-crossing rate 



27r(Do + Vo) 1/2 \l+AD b 



A 



1/2 



exp 



D b + ip b/ 
(45) 

where E (= E b — Eq) is the activation energy. In passing 
we note that the temperature due to internal thermal 
noise, the strength of the external noise, the correlation 
times and the damping constant, are all buried in the 
expression for the generalized escape rate for the open 
system through the parameters Dq, D b , ipo, ip b and A. We 
also note that k is also an implicit function of the degree 
of correlation A. 

From the structure of Eq. (|4"S")) it is difficult to under- 
stand the role of various parameters (internal or external) 
on the rate expression. We thus consider different limit- 
ing cases in the following part to see the behavior of the 
rate expression. 

First, we consider the case with no external driving 
and the internal thermal noise being <5-correlated, i.e. 

e(t) = n(t) = and (/(*)/(*') = 2r,k B TS(t - t 1 ). 



Making use of the abbreviations in Eqs.([22|) and (J24J), 
one can show that for pure (5-correlated thermal process 
<ip = ip b = 0, D b = D Q = k B T, A = ^/(jk B T), n = 
-(o_ +7), anda_ = -( 7 /2)-[( 7 2 /4)+w 2 ] 1 / 2 . Thus, the 
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general expression (|45f reduces to the classical expression 
for Kramers rate 



^0 



1/2 



exp 



E 



(46) 

Second, we consider the case with no external noise 
but the internal dynamics is non-Markovian with an 
exponentially decaying memory kernel. In this limit 
e(f) - 7r(t) = 0, (f(t)f(t')) = (V/r c )eM-\t-t'\/r c ) 
where V denotes the noise strength (V = 2 r yk B T) and 
t c refers to the correlation time of the internal noise pro- 
cesses. Then, from Eqs. ([2"2")l and (f!?5)) we obtain 

D Q = D b = k B T 

VA) = k B T (l - ^),ip b = k B T 



M = 







lb + 


(!)«-] 







~2~ T + " b 



and hence the rate becomes 



1/2 



76 

2 



exp 



E 

T7T 



(47) 

Eq. (|4"7|) is the rate expression for pure thermal pro- 
cess with exponentially decaying memory kernel and was 
derived several years earlier by Grote-Hynes [27| and 
Hanggi-Mojtabai [H|. 

Third, we consider the case where the correlation time 
of all the noise processes (internal and external) are van- 
ishingly small, 

(/(*)/(*')> = 2 7 k B TS(t-t'), 

mat 1 )) = 2D R s(t-t f ) 

where Dr is the external noise strength in the limit tr — > 
[see Eq.©] and 7 is the dissipation due to internal 
thermal noise processes. In such a case, we have 



A) = D R = k B T + — ,i> = V-> fc = 



7 



-((X- + 7) and a_ 



and hence the rate becomes 
1/2 



up 

2-KLOb 



1 

2 



exp 



1/2 



lE b 



jk B T + D B , 

(48) 

In the limit Dr — * (i.e., D e — * and D^ — > simul- 
taneously) we recover the Kramers' results (|46l) . Also, 
the degree of correlation A is implicitly present in the 
expression of Dr [see Eq.©]. In the expression (14"5|) in 



addition to T, DR/{^k B ) defines a new effective tem- 
perature due to external driving. In a different context 
where the heat bath is modulated by an external fluctu- 
ating field we have also encountered the appearance of 
the effective temperature [HI, [01 ■ 

Finally, to study the effect of degree of correlation of 
external noise on the generalized escape rate (|4"5j) . we 
consider the case where the fluctuation in the dynamics 
is only due to external source, i.e f(t) = 0. The statistical 
properties of the correlated external noises are given by 
Eq.(3). Since, in this case, the dissipation is independent 
of the fluctuations, we may assume Markovian relaxation 
so that j(t) — 7<5(t). We then obtain after some lengthy 
algebra from the general expression (|43|) 



Dr 



D, 



M = 



1+JTR+ UJ$T R 

Drtr 

1 + JTR + UJ 2 T R 

-7- (1 + 7T fl )a_, 
1 



1 + 1T R - LJ 2 T 2 ' 

Drtr 



1 



1TR 



(1 + ITr) 



2 

'2 



1/2' 



where Dr and tr are given by ([5]) and © respectively. 
Consequently, the rate becomes: 



2iroJb \ 1 + 



1 + yr R + wgrj 



JTR - U>It r 



1/2 



1/2 



x exp 



7(1 + ITr - cjgr|) 
Dr{\ + jtr) 



E 



(49) 



It is interesting to note that the expression (|49[) denotes 
the escape rate induced by pure correlated external fluc- 
tuations, which apart from depending on the strengths 
and correlation times of the noise processes, depends cru- 
cially on the degree of correlation of the external fluctu- 
ations through the parameters Dr and tr. The absence 
of thermal temperature and the appearance of the dissi- 
pation constant 7 demonstrate the non-thermal origin of 
the noise processes as well as the absence of the FDR in 
the dynamics. 



IV. NUMERICAL IMPLEMENTATION 

To judge the potentiality and applicability of our re- 
cently developed method for computation of the rate of 
barrier crossing process, we discuss here both the details 
of the working equations from the point of view of the 
numerical implementation and the corresponding simu- 
lation of our method. To check the validity and appli- 
cability of our model from the point of view of compu- 
tational implementation, we consider the dynamics in a 
bistable potential V(x) = {x A /4) — (x 2 /2) so that the ac- 
tivation energy becomes E = 1/4. We then numerically 
solve the Langevin equation |T|) by employing stochas- 
tic Heun's algorithm [29|, [3(|. The numerical rate has 
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3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 



FIG. 1: Barrier crossing rate A; as a function of dissipation 
constant 7 for two different values of effective noise strength 
Dr. The solid lines are drawn from theoretical expression, 
Eq. (|49[l and the symbols are due to numerical simulation of 
Eq.Q. The values of the parameters used are D e = D w — 2, 
r e = 7V = t = 1, A = (open circle) and A = 1 (open square). 
Inset: A; as a function of Dr for 7 = 5 and ^ A 1. Other 
parameters are same as in the main figure. 



been defined as the inverse of mean first passage time 
0, Gil HH and has been calculated by averaging over 
10,000 trajectories. In our simulation we have always 
used r e = 7v — r = 1 such that the effective correla- 
tion time tr is always equals to 1, independent of values 
of the other parameters. To ensure the stability of our 
simulation, we have used a small integration time step 
At = 0.001 so that At/r^ <C 1. To compare our theo- 
retical prediction with numerical simulation we consider 
the case where the system is solely driven by an external 
correlated colored noise [see Eq. (l4TJl) ]. As the degree of 
correlation between the two noise processes is increased 
the strength of effective noise strength Dr increases (see 
Eq.([5])). This effectively pumps more energy into the sys- 
tem through which the escape rate should increase. In 
Fig.(l) we see this effect clearly where the escape rate 
k is plotted as a function of the dissipation constant 7 
(in the limit of moderate to large friction regime) for 
two different values of effective noise strength Dr which 
are evaluated by using two extreme values of the degree 



of correlation A (0 and 1). In the inset we show more 
explicitly how the system receives more energy through 
Dr as the degree of correlation increases. This behav- 
ior suggests that in a properly designed experiment one 
can enhance the escape rate by externally controlling the 
degree of correlation between the external fluctuations. 



V. CONCLUSIONS 

In this paper, we have generalized the Kramers' the- 
ory of activated rate processes for non-equilibrium open 
system where the system is driven by two external cross- 
correlated noise processes with the assumption that the 
underlying dynamics is non-Mar kovian. The theory takes 
into account both the external and internal fluctuations 
in a unified way. The external fluctuations considered 
are stationary, Gaussian. Our treatment is valid for in- 
termediate to strong damping limit. We have shown that 
not only the motion at the barrier top is influenced by 
the cross correlation between the external fluctuations, 
it has an important role to play in establishing the sta- 
tionary state near the bottom of the source well for open 
system. The stationary distribution function in the well 
depends significantly on the degree of correlation of the 
external noise processes. We then derived the general- 
ized Kramers' rate for the open system and examined 
several limiting cases. To establish the applicability and 
potentiality of our recently developed method, we then 
numerically simulated the dynamics in a model bistable 
potential and compared with one of the limiting cases. 
Our results generated via numerical simulation reflect a 
good agreement with the corresponding values obtained 
analytically. Our numerical analysis clearly depicts that 
the escape rate can be enhanced by increasing the degree 
of correlation between the external fluctuations. 
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